Quantum optomechanics of a multimode system coupled via photothermal and 

radiation pressure force 
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We provide a full quantum description of the optomechanical system formed by a Fabry-Perot 
cavity with a movable micro-mechanical mirror whose center-of-mass and internal elastic modes are 
coupled to the driven cavity mode by both radiation pressure and photothermal force. Adopting 
a quantum Langevin description, we investigate simultaneous cooling of the micromirror elastic 
and center-of-mass modes, and also the entanglement properties of the optomechanical multipartite 
system in its steady state. 



I. INTRODUCTION 

During the last decade, cavity optomechanics has at- 
tracted the attention of a large community of physi- 
cists [H-[l] ■ Micromechanical resonators mounted on such 
systems can be cooled down to their motional ground 
state [B, 6], opening the door to the experimental study 
of quantum phenomena in mesoscopic mechanical ob- 
jects @, H| ■ Cooling is due to an optomechanical coupling 
which can be realized in several ways, exploiting the radi- 
ation pressure force @, d, t he photothermal force 1C | , 
the optical gradient force [111] . or the Doppler force 12 j 



If the coupling becomes strong enough [j] , quantum me- 
chanical correlations between the optical field and the 
mechanical resonator can be established, and robust sta- 
tionary entanglement between optical and mechanical 
modes can be generated [H, [13]. One can exploit this 
entanglement for connecting various nodes of a quantum 
network formed by an array of similar optomechanical 
systems [TBI [l6l| . 

In order to properly study this quantum regime, one 
needs a full quantum description of the optomechanical 
interaction. Differently from the radiation pressure in- 
teraction, for which such a description has been devel- 
oped time ago [l~7l . [l8j , a quantum mechanical descrip- 
tion of the photothermal force is still lacking. A semi- 
classical model of the photothermal force, based on a clas- 
sical formalism which was developed by Metzger et al. in 
Ref. [l!| , was introduced and analyzed in Ref. 20] . Sub- 
sequent works followed a similar approach to investigate 
cooling [2l|, [22j, and more recently entanglement [23j], 
in optomechanical systems where photothermal effects 
predominate. Very recently, a phenomenological model 
for describing photothermal effects associated with the 
photon-phonon-exciton interaction, has been introduced 
by Xuereb et al. 2J] , in order to explain the experimental 
results by Usami et al. that cooled a dielectric membrane 
by photothermal effects [25( . Very few system studied the 
interplay of photothermal and radiation pressure effects, 
with the notable exception of Ref. (2(| which studied the 
classical chaotic dynamics of the system in the presence 
of both effects. 

The photothermal force is generated through the fol- 



lowing steps: i) absorption of intracavity photons by the 
micro-mechanical mirror which is consequently heated; 
ii) diffusion of the produced heat through the mirror, thus 
generating a temperature gradient inside the mechanical 
resonator; iii) the temperature gradient acts as a source 
exciting thermoelastic oscillations of the internal elastic 
modes. However, the photothermal force acts simultane- 
ously with the radiation pressure force, which cannot be 
neglected in general. This latter force couples not only 
the internal elastic modes, but also the mirror center-of- 
mass (CoM) motion, to the driven cavity mode. Here we 
present a full quantum mechanical treatment of the dy- 
namics of the multipartite optomechanical system com- 
posed by the driven cavity mode, the CoM of the mirror, 
and its elastic modes, based on quantum Langevin equa- 
tions (QLE). We show that under appropriate conditions, 
the CoM and one (or more) elastic modes of the mirror 
can be simultaneously cooled down close to their ground 
state due to the coupling provided by the photothermal 
and radiation pressure interaction. For stronger coupling 
also optomechanical entanglement at the steady state of 
the system can be generated, involving either the CoM 
mode or one elastic mode. We show that optomechanical 
entanglement of the CoM mode and of an elastic mode 
with the cavity mode cannot be simultaneously enhanced 
by the photothermal force, which enhances the optome- 
chanical entanglement of the elastic mode, and tends to 
destroy the entanglement between the CoM mode and 
the cavity field. 

In the Sec. II we describe the model and derive the 
Hamiltonian of the system. Dynamics of the system is 
investigated in Sec. Ill via linearization of the QLE. In 
Sees. IV and V we exploit the steady state correlation ma- 
trix calculated in Sec. Ill to examine ground state cooling 
and entanglement properties of the system. Sec. VI is for 
concluding remarks. 



II. MODEL 

We consider an optomechanical system composed of 
a Fabry-Perot cavity driven by an intense laser, where 
one of the mirrors is a light movable mirror whose elastic 
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normal modes and CoM oscillations around its equilib- 
rium position can be treated as harmonic. The microme- 
chanical mirror absorbs the intracavity photons, and the 
resulting thermoelastic effect excites the internal elastic 
modes of the mirrors. At the same time cavity light shifts 
and deforms the surface of the micromirror via radiation 
pressure, exciting in this way both the internal elastic 
modes and the CoM of the mirror. 

In this paper, we will study the interplay between the 
photothermal and the radiation pressure force, and how 
they affect cooling and optomechanical entanglement of 
both the internal and the CoM mechanical modes. In or- 
der to properly describe the system we start from quan- 
tum elasticity theory. 

An elastic system with density p in the presence of an 
external potential V(x') is described by a Lagrangian, 
given by [27 [ 



C= dx 

IV 



ipx' 2 - V(x') - V E (x) 



(1) 



where x' = x cm + 7?.[x + u] is the position of an arbitrary 
point in the matter, in an external reference frame [see 
Fig. Q], with 1Z the rotation of the object respect to the 
reference and u the displacement vector. The elasticity 
density potential in Eq. ([T]) is given by 



Ve(x) = ^5Z«y(x)uy(x), 



(2) 



'■j 



with elasticity tensor, iiy = ^[djUi + diUj], and stress 
tensor 



( x ) =Y~P7 \- U% 3 W ~ \ 5i i Ufcfc ( x 



EySij 

3(l-2o-) 



(3) 



where Ey and a are Young modulus and Poisson ratio, 
respectively. By replacing x' into the kinetic part of the 





FIG. 1. (Color online) Scheme of the optomechanical system 
(left) and the geometry under consideration for an elastic sys- 
tem and its coordinates (right). 



Lagrangian of Eq. ([T]), one gets [26 



C =2 mx 



2 

cm 



+ \ / y ^p(x)u 2 (x) 



- / dx[V(x')+V E (x)] 
Jv 



(4) 



where m is total mass of the system and the second term 
refers to the rotational kinetic energy. In the present 
case, the system of interest is a clamped micromechan- 
ical mirror; as a consequence, the rotational degrees of 
freedom are not relevant and we will ignore them in the 
following. The Hamiltonian of the system without any 
external potential reads 



2 

Pcm 

2m 



rfx^+Vs], 

V l P 



(5) 



where the conjugate momentum of the CoM and the dis- 
placement vector are defined as p cm = <9£/9x cm and 
7Tj(x) = dC/dui(x), respectively. The equation of mo- 
tion for the deformation vector, u(x), can be easily de- 
rived from Eq. ([5]) and reads 



pu 



£yV(V.u) E Y y 2 u 

2(l + cr)(l-2cr) + 2(1 +o)' 



(6) 



One splits the vector u(x) into the transversal and lon- 
gitudinal components as u = u; + u t , where V ■ u t = 
and V x Uj = 0. Therefore, the elastic wave equations 
along the two directions read as follows, 



where 



ui - c 2 V 2 u z = 0, 
ii t - c 2 V 2 u t = 0, 



E Y (l~a) 
p{l + a)(l-2a)' 

2p(l + cr)' 



(7) 
(8) 



(9a) 
(9b) 



are the propagation speed of the longitudinal and 
transversal elastic waves, respectively. 

The photothermal effect can be included by noticing 
that when there is a temperature gradient in the material, 
the longitudinal wave equation attains a source term (27| 



U; - C 2 V 2 U; 



3p{l - 2a) 



VT, 



(10) 



where ath is the linear thermal expansion coefficient of 
the mirror. This source term can be taken into ac- 
count by introducing an effective potential, which is ulti- 
mately caused by the absorption of the intracavity pho- 
tons. This effective thermoelastic potential adds to the 
effective elastic potential associated with the trapping of 
the CoM motion, so that, neglecting for the moment the 
radiation pressure effect, the total external potential is 
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given by imwj 1 x^ m + J v ebcV p t(x'), where the photother- 
mal potential density is given by 



v ' i 

where T is temperature of the thermal reservoir of the 
mirror. The above defined V p t(x') is implicitly a function 
of external coordinates via the temperature term. Actu- 
ally, the temperature is a function of the number of the 
absorbed cavity photons, which in turn is a function of 
the cavity length (determined both by CoM and defor- 
mations on the surface of the mirror) . The distribution of 
the temperature inside the micromechanical mirror can 
be calculated from the heat equation 



K th V z T - 



ith; 



(12) 



where C is the specific heat capacity and K t h is the ther- 
mal conductivity of the material, while Qth is the heat 
flux, which in our case is due to the absorption process 
at the surface of the mirror. 

In typical situations a one-dimensional model focus- 
ing only on longitudinal deformations of the mirror 
along the cavity axis x, and neglecting the transversal 
elastic modes, provide a satisfactory description of the 
physics. The two equations above can be rewritten in 
one-dimensional form as 



Vpt(x) = 6a^) (T - To) dP 



dT d 2 T 



fit^c rout, i 

™2 ^2 



(13a) 

(x',t)8(x), (13b) 



where we have explicitly written the heat flux term, with 
I2 the intracavity photon absorption rate, /3 is quan- 
tum efficiency of the absorption, and is the radius of 
the optical mode on the surface of the micromechani- 
cal mirror, whose are is equal to S. The Dirac delta 
function assures that all the photons are absorbed at 
the surface of the mirror (23|. The explicit form of the 
effective thermolelastic potential is obtained by solving 
the heat equation of Eq. (|13b ) and inserting the solu- 
tion into Eq. (fT3"k ). The heat equation must be solved 
within the mirror, < x < i, where £ is the mirror 
thickness, with the initial condition T(x,t — —00) = To, 
and with the boundary conditions T(x = £,t) = To and 
[dT(x,t)/dx] x =o = 0. The solution is 29] 



T(x,t)=T a + 



2f3huj c 



£ cos [(2* + l)^] 



(14) 



x J* dfir (0 cx P { - [(2* + 1)^] 2 (t - f)}. 

where we have defined v 2 = -Kth/ pC '■ 



From the homogeneous elastic wave equation of Eq. ([7]) 
one can find the eigenmodes of the longitudinal elastic 
waves, forming a complete orthonormal set of functions, 
{w°(x)}, normalized so that J v dx u°(cc)u°(x) = VSij), 
where V is the volume of the micromechanical mirror. 
Therefore, an arbitrary longitudinal elastic wave, and the 
generalized elastic momentum as well, can be written as 
an expansion over this set of eigenmodes, 



u x (x,t) = 22u°(x)Q n (t), 

n 



(15) 
(16) 



where P n — mQ n . At this point, it is straightforward 
to perform a canonical quantization by promoting the 
expansion coefficients Q n and P n associated with the n- 
th longitudinal elastic mode to operators fulfilling the 
canonical commutation rules [Q n , P n ] = ih. 



A. Hamiltonian of the system 

Using Eqs. ((3J), Q13b,). and Q14p. the total Hamiltonian 
of the system can be written as 

H = H c + H m + H lp + H pt , 

where 

H c = Hu> c o/ a, (17) 
is the cavity mode Hamiltonian, 

H m = l^nA + <Zcm) + E ( P n + Qn) > (18) 

n 

is the free mechanical Hamiltonian, where we have 
rescaled all mechanical operators by making them dimen- 
sionless, such that [Q n ,Pn] = i and [<? C m,Pcm] = i and we 
have defined the elastic mode frequencies 



E y {1-(t)S 



inr 2 (l +<r)(l - 2a) J v dx 



1 du° 
dx 



9o<lc 



-J= ^2 GonQnUn(x = 0) 

V n 



(19) 



at a, (20) 



is the optomechanical interaction due to radiation pres- 
sure, and 

H pt = -fiG 0nX J2Qn I dt'h n (t-t')I2 ut (t'), (21) 

is the photothermal interaction term. In the interac- 
tion terms we have introduced the radiation pressure 
couplings go = (w c /L)^//i/mw m for the CoM, Go n = 
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{uj c / L)yj2h/mVt n for the elastic modes, with L the equi- 
librium cavity length, and the time-dependent photother- 
mal couplings 

M«)-£jf*{c£)E-[<» + i>*f] 

(22a) 



Langevin equations 



x exp | — [(2& - 

f3E Y a th LS 
3pCwr$£(l - 2a) 



(22b) 



where we have singled out the dimcnsionlcss coupling co- 
efficient x- The derived full Hamiltonian shows that ra- 
diation pressure and the photothermal force couples the 
optical cavity mode simultaneously to the mirror CoM 
and to its internal elastic modes. More precisely, the ra- 
diation pressure acts on both the CoM and the internal 
modes, while the photothermal force excites the elastic 
modes only and does not affect the CoM. 

In order to obtain an explicit expression for the pho- 
tothermal and radiation pressure couplings and parame- 
ters, one needs the explicit form of the normalized eigen- 
modes {u^(x)}, which is determined by the boundary 
conditions. In the present case both mirror surfaces are 
free to move, implying du^/dx — at x — and x = £, 
so that the normalized u^s are given by 



u° n (x) = V2 



Using Eq. (f2"2"b), one gets 



cos 1 



(23) 



M*) = E 



k=0 



4n 2 exp{-[(2fc + l)ff] 2 t} 
4n 2 -(2fc+l) 2 



(24) 



showing that in practice, the characteristic time of the 
photothermal force is determined by the largest thermal 
diffusion time, that is r t h = (2(,/nv) 2 . Moreover, from 
Eq. (|I9p we arrive at the following explicit expression for 
the elastic resonance frequencies 



fif. 



E Y (l - a)S 



pnr z {l +cr)(l - 2a) 



nir\- 
~T) 



(25) 



III. DYNAMICS OF THE SYSTEM 

Dynamics of the system can be fully characterized by 
its QLE, which can be derived from the Hamiltonian 
derived in the previous section and including damping, 
noise terms, and also driving of the cavity mode through 
the input mirror with coupling rate K\, by a laser with 
frequency ojl and power V . In the frame rotating at the 
laser frequency one has the following set of 27V + 3 of 



Pcm 
Qn 
Pn 



-^m<7cm ~ 7mPcm + 9vo) CL + f , 



£1 P 



—£l n Q n — Y n P n + Go n a) a 

+G nX f dt'g n (t - i'K 2 ° ut (0 + S n , 



(26a) 
(26b) 
(26c) 

(26d) 
)a + S 
(26e) 

where in writing Eqs. (J26jl) and (f26b) we have exploited 
Eq. dHJ} to get u° n {x = 0) = y/2. We have also defined 
Ao = lu c — cjl and £ = ^/2k{P ' /Ti/jj-^. The noise operators 
of the system are zero-mean Gaussian noises; a™ and aif 
are the vacuum input noises of the cavity field entering 
from the two mirrors, one with decay rate K\ and the 
other with k 2j so that the total cavity decay rate is k c = 
K± + K2- The correlation function of the two noises is 
given by 



-(« c + iA )a + i(goq cm + ^ G 0n Qr 

n 

-v / 2KTa i 1 n + ^/2K2"4 n , 



(af(ty»'\t')) = [N(lo c ) + l]6i S S{t - t% 
(a}^(t)af(t / ))=N(u Jc )S lJ S(t~t / ), 



(27a) 
(27b) 



where i, j = 1,2 and JV(w c ) = ( exp{fkd c /kQT} — l) is 
the mean equilibrium thermal photon number (&b is the 
Boltzmann constant and T is temperature of the reser- 
voir). At optical frequencies hw c /k-QT 3> 1 and therefore 
N{u>c) ~ 0, so the only relevant correlation function is 
that of Eq. (|2~Ta) . The mechanical modes (CoM and the 
internal elastic modes) are affected by a friction force 
with damping rates 7 m for the CoM mode and T n for 
the elastic modes. These viscous forces are associated 
with Brownian stochastic forces obeying the correlation 
function [301 



mm) 



Jul I ^ e -*-(*-f) w 



2vr 



i -<■<■: I. (rr^^a) 



-2/«p 



(E k (t)E k ,(t'))=5 kk ^ f p 
ilk J 2n 



e -iu(t-tf) u 



1 + coth ( 



•2fc B T^ 



(28b) 



where k,k' = 1, N. 

We are interested in strong optomechanical coupling, 
which is the optimal regime both for cooling and entan- 
glement. However, the single photon couplings go and 
Go n are typically of the order of few KHz, and there- 
fore they are not able to provide the required coupling 
strength in the weak field regime. One can circumvent 
this problem by increasing the intracavity field amplitude 
(a), obtaining in this way an effective linearized large 
coupling strength. By increasing the input power and 
employing high-finesse optical cavities, high intracavity 
intensities can be achieved. At this situation, the cavity 
field is described by an intense coherent state and every 
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mechanical mode is shifted to a new equilibrium posi- 
tion depending upon the number of intracavity photons. 
Therefore, every operator of the system can be decom- 
posed into two parts: a classical steady state value, and a 
fluctuating quantum part around this steady state value, 
that is, for a generic operator O we have O = (O) + SO. 
By replacing in the nonlinear QLE of Eqs. ([2"B)) . we get 
the steady state values as the following 



(?cm) = 3o(a) 2 Mi 

(Qn) = ( 1 



where we have introduced the effective cavity detuning 
A c = A - go(q C m) - J2n G 0n(Qn), and we have cho- 
sen the phase of £ can be tuned so that (a) is real and 
positive. 



where we have introduced the cavity field quadratures 
SX a and 5Y a , so that we have 5a = (SX a + iSY a )/y/2 and 
their corresponding noise operators from each side of the 
cavity X£, X% Y£, and Y$. The effective optome- 
chanical couplings between the quantum fluctuations are 
given by g = go(a)\^2 and G n = Go n (a)^/2. Because of 
the integral appearing in Eq. (JUJi), we need to solve these 
equations in frequency space using the convolution the- 
orem. In the frequency domain, we can write the above 
equations in the following compact form 



- \iuil + A(u>)] u = n, 



(30) 



where I is the (2iV + 4) x (2iV + 4) identity matrix and we 
have defined respectively the vector of system operators 
u and the noise vector n as 



A. Dynamics of the quantum fluctuations 

For the small quantum fluctuations, we arrive at the 
following linear equations 



Sq, 
Sp, 



cm 
'cm 



J cmi 
5cm 



Ym^Pcm + gSX a + £, 



SPn = 



sx a 

5Y„ 



ui m Sp, 
-cj m Sq ( 



(29a) 
(29b) 
(29c) 



sx a 



+ X I dt'h n (t - t')[2K 2 8X a - V2^X^9d) 

/2^X329e) 

^y£ 

(29f) 



-n c SX a + A c 5Y a + V2^K n i 
- Kc SY a - A c 6X a + V2i^Y^ - 
+gSq cm + G n SQ n , 



J 



u = [Sq cm ,Sp cnu 6X a ,5Y a ,6Q 1 ,SP 1 ,SQ 2 ,SP2, ■ ■ -] T , 

(31) 

n = [0, 1 V2^X™ + s/2^X? a , V^iYl: + sf& 2 ™, 

0,-V2^xGih 1 (oj)X i 2 l + E 1 ,---] T , (32) 

where h n (uj) is Fourier transform of the response func- 
tion. By keeping only the term associated with the 
largest thermal diffusion time, it is straightforward to 
get 



h n (u>) 



4n 2 r th 



A(co) 



(l-4n 2 )(l-^T th )' 
The drift matrix A(u) is given by 



(33) 



( 


w m 






















-Win 


~7m 
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A c 














g 







-Ac 




G x 





G 2 






















Qx 














Gi(l- 







-Qx 


-Tx 































fi 2 








G 2 (l- 


1- 2k 2 X^2(w)) 











-n 2 


-r 2 



(34) 



The steady state associated with Eq. (f30|) is reached 
when the system is stable, which occurs if and only if 
all the eigenvalues of A(lu — 0) have negative real part. 
These stability conditions can be obtained, for example, 
by using the Routh-Hurwitz criteria. In this paper, we 
shall restrict to the situation with A c > 0, i.e., with a 



red-detuned laser, and in this parameter region the only 
nontrivial stability condition is 

A c r g 2 Gl / 8n 2 \1 „ 

(35) 
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B. Correlation matrix 

Since we have linearized the equations around the 
steady state values and all the noise operators of the sys- 
tem are Gaussian with zero-mean, the steady state of 
the quantum fluctuations is a zero-mean Gaussian state, 
whose properties are fully characterized by its correlation 
matrix. This steady state correlation matrix is defined 
as 

Vij = -(ui(oo)uj(oo) + ttj(oo)ttj(oo)). (36) 



from the formal solution of Eq. (|30|) one gets the following 
expression for the steady state correlation matrix 



V 



duM(uj)D(uj)M^{uj) 



(38) 



Since at a generic time t one can write 



KM) 



4-7T 



,-i(ui+ui')t/~ 



(37) 



where M{uS) = [iul + A(w)] -1 , and we have de- 
fined the diffusion matrix D{uj) through the relation 
(rii(uj)fij(uj') + n J -(u; / )n i (u;))/2 = D(w)8{u + uj'), whose 
explicit expression is 






























7 m u; r hjj -, 
coth — 






















uj m L 2fc B T J 




























-K 2 Gixhi{u)) 
















Kc 












































-K 2 Gixhi(uj) 








CiM 


















































CaH 



and where Cn(w) an d y>ij(ui) are defined as 

Cn(w) = K 2 (G nX \K{u:)\) 2 + ^wcoth [J^L] , (40) 
Vii (w) = 2 K2X 2 G J G,3fi[/i l (w)/ l *( w )]. (41) 

IV. SINGLE ELASTIC MODE 



A. Cooling 

Ref. [13] has shown that photothermal effects can im- 
prove cooling of a single internal elastic mode, compared 
to the situation in the presence of radiation pressure only. 
However, radiation pressure always couples the cavity 
field also to the CoM mode, and therefore for a com- 
plete description of cooling one has to include also the 
CoM, and to consider the simultaneous presence of the 
two kinds of mechanical modes. 

In Fig.[2]the effective phonon number of the CoM mode 
(blue full line) and of the elastic modes of interest (red 
dashed line) is plotted versus the effective cavity detun- 



In order to see the effects of the interplay between the 
radiation pressure and photothermal force, we consider 
a situation similar to that of Ref. (23|, i. where the de- 
tection bandwidth is chosen so that it includes the CoM 
mode and one elastic mode only. In this section we focus 
in particular on the simultaneous cooling of the elastic 
mode of interest and of the CoM mode, and also on the 
entanglement properties of the tripartite optomechanical 
system. We choose an achievable set of system parame- 
ters, which are listed in Table |TJ 



TABLE I. Parameters of the optomechanical system. 



Parameter Symbol Value 

CoM frequency cj m /27r 20 MHz 

mechanical quality factor Q m 10 5 

mirror mass m 5 ng 

cavity length L 1 mm 

laser wavelength Al 810 nm 

reservoir temperature T 4 mK 
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FIG. 2. (Color online) Effective phonon number versus cavity 
detuning in the good cavity limit n\ — 0.05o; m (left panels) 
and in the bad cavity limit K\ = 5w m (right panels): (a) 
K2 = 0, (b) At2 = 0.5/ti, and (c) k% = n\. The elastic mode 
frequency is fii ~ 0.54cj m , r th = fi^ 1 ,and \ — 0.13. Power 
of the laser is V — 15 mW. 



ing A c , for two different input rates K\, corresponding 
to two completely different regimes: the resolved side- 
band regime {n\ = 0.05w m , left panel) and the bad cav- 
ity regime (ki = 5w m , right panel). We consider three 
photothermal strengths: (a) K2 = (only radiation pres- 
sure), (b) K2 = 0.5ki, and (c) n-i = K\. From the hgure, 
one can easily conclude that including photothermal ef- 
fects considerably influences the cooling of the mechani- 
cal modes of the system both in the good and in the bad 
cavity regimes. In fact, while slightly improving cool- 
ing of the elastic mode, photothermal force worsens the 
ground state cooling of the CoM. This may be explained 
in terms of the increase of the cavity mode bandwidth 
caused by photon absorption in the micromirror. 

This effect is significant in the bad cavity limit, so that 
we cannot even see the curve of the CoM phonon number 
in right panel of Fig.^c). In Fig. [21 we have also plotted 
n e g of the elastic mode in the absence of CoM mode 
(same as the case of Ref. [HI). In the case of n c ^> uj m 
the absence or presence of the CoM mode does not change 
the situation, so the two curves coincide. However, in the 
resolved sideband limit, one can see that the presence of 
the CoM mode slightly improves cooling of the elastic 
mode toward its ground state. 

In Fig. |31 we present the density plot of the effective 
phonon number of the CoM mode and of the elastic mode 



versus laser input power and absorption rate when the 
cavity detuning is set to A c = u) m (top panels) and 
A c = £li (low panels). For the CoM mode, increasing 
the laser power improves the cooling process; at fixed 
power the way photothermal force affects cooling of the 
CoM mode depends upon the cavity detuning: for a cav- 
ity tuned to w m the bolometric force has a detrimental 
effect on cooling (Fig. G^a)), while, for a cavity tuned to 
fii, that is to the motional sideband associated with the 
internal mode, cooling improves for increasing absorption 
losses. However the relevant parameter for simultaneous 
cooling is just the total cavity decay rate n c . In fact, the 
interval for the detuning A c within which one has a sig- 
nificantly low value of n e g for the CoM mode is given by 
— k c < A c < w m + k c , consistent with the expression 
of the net laser cooling rate in the presence of radiation 
pressure interaction [3ll433| |. By increasing the absorp- 
tion rate, the total cavity decay rate n c increases, and 
the mechanical modes falling within this bandwidth are 
cooled, even for small laser powers (see also Ref. [HI). 

In Fig. |4] we show n c g for the CoM mode and for a 
selected internal elastic mode as a function of the inter- 
nal mode index n, and therefore of the elastic resonance 
frequencies f2i(n). Blue circles refer to the case when 
the optomechanical interaction is provided only by radi- 
ation pressure, while green crosses and red squares re- 
fer to two different photothermal strengths, K2 — 0.5ki 
and k% = Ki, respectively. On the one hand, cooling for 
the CoM mode, is achieved both with and without the 




00 05 0.10 0.15 0.20 0.25 0.00 0.05 0.10 0.15 0.20 0.25 








'0.1 










/A 



0.00 0.05 0.10 0.15 0.20 0.25 



0.00 0.05 0.10 0.15 0.20 0.25 
K 2 /0J m 



FIG. 3. (Color online) Density plots of n e ff as a function of «2 
and laser power, (a) and (c) refer to the CoM mode, (b) and 
(d) to the elastic mode. The upper panels refer to A c = u m , 
while the lower panels are for A c = Q 1 . The cavity input 
rate is «i = 0.05a; m . The other parameters are the same as 
in Fig. H 
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FIG. 4. (Color online) Effective phonon number of the CoM 
mode (a) and of the elastic mode (b), as a function of the 
elastic mode frequency for three different photon absorption 
rate values: K2 = (blue circle), k,2 = 0.5ki (green cross), 
and K2 = (red square). The detuning is set to A c = oj m 
and x — 0.03. Laser power is V = 15 mW and k\ — 0.05tj m . 



photothermal effect, provided that the elastic mode fre- 
quency is not to close to the CoM frequency u> m . On 
the other hand, cooling of the elastic mode improves for 
increasing absorption losses K2 7^ 0, even though is op- 
timal only for elastic resonance frequencies not too far 
(and not too close) from the CoM frequency uj m . Within 
there two narrow intervals around w m in fact, the elastic 
mode frequency satisfies the optimal resonance condition 
A c ~ fii(n), but it is not too close to uj m where a de- 
structive interference phenomenon between the two cool- 
ing processes takes place preventing simultaneous cool- 
ing, as explained in detail in [34j . In the present case 
this occurs when f2i(n = 37) ~ w m , where cooling is 
practically absent. 



B. Entanglement 

We now study the entanglement properties of the tri- 
partite system formed by the CoM mode, a selected elas- 
tic mode, and the cavity mode, by focusing on the en- 
tanglement of the three possible bipartite systems at the 
steady state. As a measure of bipartite entanglement 
we choose the logarithmic negativity [3oT - l37| . given by 
En = max[0, — In 2^"] , where rf~ is the minimum sym- 
plectic eigenvalue of the partially transposed bipartite 
covariance matrix (3o| . 

We first investigate the dependence of the logarithmic 
negativities upon the effective cavity detuning at differ- 
ent values of the cavity absorption losses K2 (see Fig. [SJ . 
In Fig. [5] we have also included for comparison the neg- 
ativity of the elastic mode-cavity mode bipartite system 
in the absence of the CoM mode (black dot-dashed line) . 
Since the stability condition in this latter case (i) dif- 
fers from the stability conditions of the tripartite system 
including also the CoM mode (ii), we have used differ- 
ent colors — dark purple for case (i) and light purple for 
(ii) — to denote the unstable region in each case. One can 
see from Fig. [5] that when there is no photon absorption 



(K2 = 0), the presence of CoM mode does not appreciably 
affect elastic mode-cavity mode entanglement, and only 
the instability of the system prevents this bipartite op- 
tomechanical entanglement to achieve values large than 
those achieved without the CoM (i). In the presence of 
a nonzero photon absorption rate, decoherence of both 
systems is increased, and therefore the bipartite -Bat's 
decrease (see Figs. 5(b) and (c)). However, while for the 
case without the CoM (i) the stable region rapidly ex- 
tends, allowing the Em of the elastic mode-cavity mode 
subsystem to reach higher values, in the presence of the 
CoM (ii) the stability region is not significantly altered 
by the photothermal force, and the En cannot reach sim- 
ilar values (see Fig. EJb)). At the same time a small but 
nonzero CoM-cavity mode entanglement in a narrow in- 
terval of detuning is found. Further increase of photon 
absorption K2 has a detrimental effect of the bipartite 
optomechanical entanglement, so that there remains no 
CoM-cavity mode entanglement when k% = K\, while the 
elastic mode-cavity mode entanglement decreases even 
though it remains non-negligible. We notice in addition 
that the CoM mode-cavity mode En, differently from 
the elastic mode — cavity mode En, is not maximum at 
the bistability threshold, and that by increasing the pho- 
tothermal force the first region at which the entanglement 
disappears is just at the edge of instability region. 

We finally notice that for the parameter region con- 
sidered here, there is no entanglement between the two 
mechanical modes. Here we found results consistent with 
those of Ref. [34| , which showed that small but nonzero 
bipartite mechanical entanglement can be obtained only 
in a regime where the oscillators are heavily damped and 
the cavity finesse is very high (k c ~ 7 m < u m ). 

To see that to what extent the photothermal force in- 
fluences bipartite optomechanical entanglement and to 
examine its relation with the stability conditions, we plot 
En of the CoM mode-cavity mode and elastic mode- 
cavity mode bipartite subsystems as a function of the 
absorption rate K2 and of the power of the pump laser 
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FIG. 5. (Color online) Logarithmic negativity of the three 
distinct bipartite subsystems: CoM mode-cavity mode (blue 
solid line), elastic mode-cavity mode (red dashed line), CoM 
mode-elastic mode (green dotted line), and elastic mode- 
cavity mode when the CoM mode is absent (black dot-dashed 
line), as a function of A c ; (a) K2 = 0, (b) K2 = 0.5ki, and 
(c) K2 = Ki- Laser power is V — 50 mW, Ki = 0.25u m , 
fii ~ 0.54xj m , x — 0.13, and r t h = 1/fii. 
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FIG. 6. (Color online) Density plot for the logarithmic 
negativity of CoM mode-cavity mode (a) and (c) , and of 
elastic mode-cavity mode bipartite subsystems (b) and (d), 
as a function of k,2 and laser power. Cavity detuning is 
A c = 0.4o; m for the upper panels and A c = 0.9cj m for the 
lower panels, and the cavity input rate is ki = 0.05o; m . The 
other parameters are the same as in Fig. [5] 



[En of the bipartite subsystem composed of the mechan- 
ical modes is again zero in the whole parameter region). 
From Fig. EJa) and (b), which refers to A c = 0.4w m , 
one can easily see that the presence of photothermal ef- 
fects leads to a narrowing of the instability region. As a 
consequence, while the maximum value of elastic mode- 
cavity mode Em is about 0.21, when K2 = 0, by increas- 
ing the photon absorption rate, the instability threshold 
is "pushed back" and the elastic mode-cavity modcEM 
may reach values as high as 0.26. On the contrary, the 
CoM mode-cavity mode Em at this value of the detun- 
ing A c is always absent. Instead, for A c = 0.9u; m which 
refers to Fig. |^c) and (d), one can achieve a small value 
of the optomechanical entanglement for both the CoM 
and the elastic mode, at high pump powers. We notice 
from Fig. Etc) and (d) that there is an optimal nonzero 
value of K2 for both the CoM mode-cavity mode and elas- 
tic mode-cavity mode entanglement giving a maximum 
value for Em- This fact occurs also for the CoM mode 
even though it is not affected by photothermal effects, 
because the entanglement depend upon the total cavity 
decay rate k c = k± + K2 which must be not too small in 
order to have an appreciable entanglement. 

Similarly to what we have done for cooling, we fi- 
nally analyze the behavior of bipartite optomechanical 
entanglement as a function of the difference between 
the CoM and elastic mechanical frequencies. In Fig. [7] 
we plot CoM mode-cavity mode (Fig. [Tfa)) and elas- 
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FIG. 7. (Color online) Logarithmic negativity of CoM mode- 
cavity mode (a), and elastic mode-cavity mode (b) bipartite 
subsystems as a function of the elastic mode frequency. The 
detuning is set to A c = 0.87aj m and \ — 0.03. Laser power is 
V = 100 mW and m = 0.25o; m . 



tic mode-cavity mode (Fig. [T^b)) entanglement versus 
the elastic mode index number which determines the 
value of ft i, for three different values of the intracav- 
ity photon absorption rate K2- Similarly to the re- 
sults of Ref. 



34J, when the frequency of the elastic 



mode is close enough to the CoM frequency (f2i(n 
37) ~ Wm), the two optomechanical entanglements van- 
ish. In fact, as illustrated in Ref. (34j, when the two 
mechanical modes are at resonance, the cavity mode is 
strongly coupled, and entangled, to the collective coordi- 
nate (go?cm + GoiQi)/(<7o + G'oi), and uncoupled from the 
relative coordinate (goq C m~GoiQi) / (go+Goi) ■ The CoM 
mode and the single elastic mode are linear combinations 
of these two coordinates and therefore their entanglement 
depends upon the state of these coordinates. At T = 0, 
the CoM mode and the elastic mode of interest are en- 
tangled with the cavity mode, due to the strong coupling 
of the collective coordinate with the cavity mode and be- 
cause the relative coordinate, even though uncoupled, is 
in its ground state and does not degrade the established 
optomechanical entanglement. However, as soon as tem- 
perature is increased, the thermal fluctuations of the un- 
coupled relative coordinate prevails over the effect of the 
optomechanical coupling and destroys the entanglement 
of the CoM and of the elastic mode with the cavity field. 

Furthermore, it is interesting to note that for the 
elastic modes at optimal distance from the CoM mode 
(31 < n < 34 and 40 < n < 43) a the presence of 
photothermal coupling gives higher optomechanical en- 
tanglement with respect to the situation with only the 
radiation pressure interaction. Fig. [7] also show that en- 
tanglement sharing also occurs in the steady state of 
the tripartite system under study. In fact, the elastic 
mode — cavity mode entanglement is redistributed to the 
CoM mode — cavity mode entanglement as soon as radia- 
tion pressure coupling and therefore the CoM is present. 
More in detail, when Cl± < w m the elastic mode-cavity 
mode Em is much greater than CoM mode-cavity mode 
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entanglement, while for fix > oj m the opposite situation 
occurs. Hence, one can increase the optomechanical en- 
tanglement of one mode at the expense of the other by 
selecting the appropriate internal elastic mode. 



V. CONCLUSION 

In conclusion, we have introduced a general Hamilto- 
nian including both photothermal and radiation pressure 
interactions in an optomechanical system formed by a 
Fabry-Perot cavity with a movable micromirror. Even if 
referred to this specific system, one could apply this anal- 
ysis to different optomechanical devices and geometries. 
The two kind of interactions couple a driven cavity mode 
to the CoM mode of the micromirror and to several in- 
ternal elastic modes, and we have analyzed the quantum 
properties of this optomechanical system, by adopting a 
QLE approach. By linearizing the QLE of the system, 
cooling of the mechanical modes has been studied: we 



have shown that in addition to the CoM mode, one can 
cool the internal elastic modes of the mirror toward their 
ground state. This simultaneous cooling is helped by the 
photothermal force in the sense that the parameter re- 
gion where simultaneous cooling is achieved is extended 
by its presence. As first shown in Ref. [34j . cooling van- 
ishes because of a classical destructive interference only 
when the CoM and the elastic mode mechanical frequen- 
cies becomes close to each other. We have then studied 
bipartite entanglement of the tripartite system formed by 
the CoM mode, a selected internal elastic mode and the 
cavity mode, focusing especially on the two optomechan- 
ical bipartite systems, because entanglement between the 
two mechanical modes is achievable only in a very narrow 
set of parameter space. As expected due to entanglement 
monogamy, the optomechanical entanglement of CoM 
mode and elastic modes with the cavity mode cannot 
be simultaneously enhanced by the photothermal force. 
In fact, the bolometric force enhances the optomechan- 
ical entanglement of the elastic mode, while destroying 
the entanglement between the CoM mode and the cavity 
mode. 
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